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1. INTRODUCTION 


The intensive study of elect rom a;.;n cl i c wave propagation through turbulent mediums 
began in Iho early ).!)50's. The fi '•s', work of .significance was performed in Russia 
by Kolmogorov . Chernov", and later Tntarski* . Although ti great deal of work in this 
area was performed in the Uni'cd States, intensive study did not begin until the early 
lStGO's. One of the major driving forces which spurred interest In wave propagation 
through a turbulent medium wan the advent of the laser, with its pi omise of data trans- 
mission bondwidths many ord r rs o; magnitude higher than possible with the highest 
f req u e nc y i n i c r own vo s . 


Al' hough initially there - as some interest in high-datn-rnto trt is mission between points 
on the cart! 's ■ urfaco, it soon became apparent that the effects of weather and the 
line-of-sight nature cf laser transmission would preclude high-data-rate laser com- 
munications except over short distances and under favorable weather conditions. The 
situation, however, for ground-to-spacc and space-to -ground communications is much 
more favorable, and therefore, work on this type of data link has been actively pursued 
by both NASA and the Air Force. Studies* have shown that by the use of redundant 
earth-based transmitter/ receiver sites, link outages can be reduced to a few tenths 
o r a percent. Through che use- of airborne transmitter/ reeel vor systems, even these 
minor outages can be eliminated. 

From the beginning of the study of atmospheric turbulence, the sophistication of Do 
theoretician has far outpaced the experimentalist. There are good and sufficient 
reasons for thin state of affairs. The mechanism for the generation of turbulence is 
beyond the control of the experimen'er and, therefore, it is impossible to perform con- 
trolled oxportmmts in which the state cf turbulence can be fixed while measurements 
are performed. The problem is further exacerbated by the fact that the experimenter 
rarely has sufficient instrumentation to even specify the state of the atmosphere during 
his measurements. 


Some of the earliest experiments in optical propagation wore made by astronomers, 
who wore concerted with the variations in seeing and scintillation because these effects 
masked their abil.ty to study the characteristics of extraterrestrial sources. Although 
considerable timo and effort wus expended in their studies, most of their work was 
performed during a peri;>d when no adequate theory was available to explain the mech- 
anism of the obsetved effects, and therefore much of the work Is of n subjective nature 
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with inadequate supportive data to interpret tho results. Notable exceptions :o this 

G V B 

state of affairs arc the work of Mikescll and Hrag , Protheroe , Young , Kaniscy , 

9 10 

Coulman , and later, Pufton . Except for the studies made using astronomical 
sources, nearly all of the experimental attempts to verify theory have been made 
over horizontal propagation ranges. If a review of the proportion of work performed 
on horizontal ranges compared to vertical propagation paths wore made, it would re- 
veal that while thousands of experiments have been performed over horizontal ranges, 
almost nothing has been done to measure the effects of atmospheric turbulence on a 
laser beam pronngating ove: a vertical path. 

When work on the study of atmospheric turbulence was started at GSFC, by tho nutuors, 
it became apparent after a short amount of study that the theoretical predictions of 
turbulence effects were being adequately pursued by other agencies and organizations 
and hat the area which was most in need of investigation was the experimental verili- 
cation of existing theory. Because all of the theoretical predictions of tho optical 
effects are based upon th«> slate of ntmr. :pherlc turbulence, n development procram 
was started to build a sensory device capable of directly measuring the state of tur- 
bulence. The outcome of this effort was the thermosonde, a lightweight electronic 
device capable of being borne aloft by a small weather balloon to measure the state of 
turbulence at altitudes to 100,000 feet, well above the lust significant turbulence. In 
1908, interest in laser propagation was at it., peak, and many groups were considering 
methods of measuring turbulence in the "pper atmosphere, ltceniise of their ability to 
reach high altitudes *’nd because they generate no atmospheric turbulence, the authors 
were at this lime considering mensyjarffumts mode from high-altitude research balloons. 
Because of our interest in balloons os a vehicle for turbulence measurements, the 
Oflico of Advanced Research and Technology forwarded to us a proposal by the Sylvanin 
Electronic Systems, Western Division to perforin n series of balloon-borne laser 
propagation experiments. The outcome of this proposal was the experiment described 
In this report. 


The experiment, which was Titled "Bnlloor Atmospheric Propagation Experiment" 
(BAPE), was performed at White Sands Missile Range In the fall of 1970, and is to 
this date the only controlled measurement of laser propagation over a vertical path. 

The experiment was repented in 1971 at tho same location and the two flight series havo 
been named BAPE I and BAPE II to distinguish between them. This report will descrlbo 
only the data obtained in BAPE 1. 
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The experiment plan was to measure the optical effects of turbulence over a vertical 
oath while oi the same time measuring, with the thermosonde, the distribution of 
turbulence along the path. In this manner, the optical effects predicted by theory 
could be compared to the measured restliis in an effort to verify the validity of 
theory. The experiment has not uncovered any phenomena in contradiction to existing 
theory, but has confirmed many of the existing predictto'.s. We feel that although some 
new and unexpected phenomenon might perhaps be n.ore dramatic, our confirmation of 
existing predictions is actually more constructive. 

To this date there nave been only t.rcc successful attempts to measure the character- 
istics of laser propagation between ground and space. All of these hnve been per- 
formed by the authors, starting with the GEOS II Laser Detector Experiment and fol- 
lowed by the BA PE 1 and BAPE 11 experiments, .he results of these vxperiments con- 
firm me predictions of existing theory withi:. the limits of the experiment Because of 
the importance of this field of investigation to the development of the laser as a corn- 
man. cations tool and because of the demonstrated success of balloon-borne instrumen- 
tation as a means of obtaining atmospheric propagation data, future experiments of 
ti..s type ore being planned. 

The balloon-borne experiment described in the subsequent sections of this report con- 
stitutes the lirst definitive set of laser data on iho effects of atmospheric turbulence 
over a vertical or near-vertical path. As such, the experiment provides the basis lor 
meaningful predictions of laser system performance and data which establishes the 
utility and limitations of various propagation theories. Therefore, this was an impor- 
tant experiment with far-reaching consequences. 

Basically the objective of the experiment was to propagato two laser beams, one argon 
(bl-s. r >A°) and one CO, (10. Gp), over a vertical path from the ground to receivers lo- 
cated above the atmosphere, and to measure the scintillation of the upgoing beams. 
Simultaneous measurements oi the slate of turbulence along the path allow the scin- 
Ullatkm statistics to be compared on the basis of existing theory to the experiment 
results. This type of experiment has four important results: 

1. The measurements of scintillation arc directly applicable to estimation of 
scintillation of earth-to-spacc laser links. 

2. The measurements of the stale of turbulence protide now data oo the state 
of atmospheric turbulence, no direct experimental data on turbulence above 
2.5 1 m (10,000 ft) was available prior to the experiment. 
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By comparing the nclunl scintillation .statistics with those based upon 
theory plus the turbulence profile, we can check the validity of existing 
propagation theories . 1 

Knowing that these theories are correct (assuming that the theories check 
out), ue can make far ranging predictions on the performance of laser 
eo • inimical ions systems with a high degree of confidence. 
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2. SUMMARY 


This section provides n brief synopsis of the results from the Balloon Atmosplx ric 
Propogntion Experiment \ 

The vertical profile of turbulence in the first C km of the atm ■'sphere was found to be 
approximately 

„2 . ,.-14. -1/3 -h/1000 -2/3 

Cu = 4 x 10 h c m 

Thelog-amplitude variance at 45° from zenith had an average value of 0. 1!)8 tit 
0.5145 pm and 0.00510 at 10.0 pni. Measured wavelength dependence of scintillation 
was found to be within 8'.c of that predicted by the 7/G power law. The sample log- 
amplitude variance was itself a random variable with a log-normal distribution and a 
standard deviation of ln2. This indicates that the strength of turbulence varies with 
time in the same manner. 

Log-amplitude covariance was measured and the correlation lengths were found to be, 
as predicted by theory, proportional to K ' the s-tc ol the c irrelntion lengths was 
also found to be as predicted by lheo*v. In order to test this last relationship, a 
method of approximating the atmosphere by a multilayer system of turbulence was 
developed and used to predict the correlation lengths. Both correlation length and 
the shape of the correlation coefficient predictc ’ by this procedure show close cor- 
respondence to measured values. 

Using the above method, aperture averaging for large apertures is computed and shown 
to be in close agreement with other methods of obtaining aperture averaging factors. 

Autocorrelation of irradiancc and power spectral density were measured and compared 
to predicted results. A close agreement between theoretical and measured values indi- 
cates that t he tempor al behavior o r the scintillation can be obtained from a knowledge 
of the height of the turbulence layers and their wind speeds. 

The results of the experiment represent a fairly complete investigation of the scintil- 
lation effects of atmospheric turbulence over a vertical path, although much remains to 
be measured. We have found no parameter which does not have the magnitude and 
t : mvior predicted by theory. Results of the experiment tend to confirm the validity of 
present predictions of the effects of atmospheric turbulence. 
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3. DESCRIPTION OF EQUIPMENT 


3.1 RAPE PAY I. PAD 

The following sections describe the design anti construction of the RAPE payload which 
was the flight vehicle used to carry the optical detectors of the experiment above the 
atmosphere. The payload was fabricated bv Sylvar.la Electronic Systems Western ■ 
Division under contract NAS 5-jl»i7‘J. The material from the following sections is in 
part abstracted from the RAPE Phase I report of August 14, 19(511, with changes and 
additions made where required. The principal authors of this report w>-re Mr. Arthur 
Kraemor, Dr. Paul J. Titterton, Mr. Tom V.ulcahy , and Mr. beo.t Overstreet. Other 
contributors include Mr. Clyde Drown, Mr. Connell Ward, Mr. Miel.-.el Pnlus, and 
Dr. Douglas Woodman. Professor J. Richard Kerr of the Oregon Graduate Center 
acted, as a program eons'.!'. 1 . a n * . 

The pm'|>oso of the material in the following sections is to familiarize the reader with 
the construction and operation of the payload so that he may belter understand how the 
data was obtained and its significance. 

3.1.1 GENERAL DESCRIPTION OF PAYLOAD 

The RAPE payload shown in Figures 3-1 and 3-2 consisted of an optical receiver sys- 
tem W an inverted azimuth-elevation gimbal system. The elevation axis (see Figure 
3-3) contains four optical receivers, a star tracker, servo electronics, and two cube 
corner retroreflectors. It is pivoted on bearings supported by the azimuth gimbal and 
driven by a dc torque r through elevation angles from -2(5 to -90 degrees. The azimuth 
gimbal supports the elevation axis and connects it to the azimuth axis. The large box- 
shaped section of the azimuth gimbal contains batteries, telemetry encoders, raid other 
electronics for the operation of the servo system. These electronics arc connected by 
a cable wrap to the elevation axis. Above the azimuth gimbal is the azimuth bearing 
and slip ring assembly which provides continuous rotation in the azimuth direction and 
electrical connection to the gondola. A dc torquer, tachometer, clutch, and azimuth 
encoder are attached to the top of the azimuth bearing and slip ring assembly to com- 
plete the azimuth drive system. The gondola (see paragraph 3. 1 -5) provides a frame- 
work for supporting the components attached below the azimuth axis and for supporting 
the battery boxes. It consists of three arms connected at their inner ends to the ele- 
vation axis and at their outer ends to the battery boxes. The battery boxes contain the 
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i, Figure 3-2. BAPE Payload 







batteries necessary to power the payload as well ns telemetry and electronics systems. 
The mass of the heavy batteries on liit> ends of the arms acts as pn inertia flywheel 
against which the n?imuth axis drives. Struts connected to the bottom of the battery 
boxes provide protection of the payload when it strikes the ground at the end of the 
parachute descent. Disposable, cor regaled cardboard crush pads on tr.o bottom of 
these struts are used to absorb the energy of the impact. 

3.1.2 ELEVATION AXIS AND OPTICAL DETECTOU SYSTEMS 

An optical sehemaiic of the elevntion axis and optical detector systems is shown in 
Figure 3-4, The elevation axis contained two detectors at each of the wavelengths 
0.51:3 pm and 10.0 pm. Due system o|>o rated fixed and the second system operate 1 
with a movable input mirror. Iti system I, the received beam was received by the di- 
chrolo beam splitter, which was used to separate the argon (0.51-13 pm) and the CO., 
(10.0 pm) beams. The Uchroic reflected the 10.0-pm beam through a refractive f/1 
germaniui • lens onto the IlgCdTe detector mounted in a liquid nitrogen dowar. The 
0.5145-pm beam passed through the dlehrole and was reflected through a lens arrange- 
ment onto the photomultiplier detector. 

In system II, 'he received beam was reflected off the movable input mirror and a . ixed 
mirror prior to reaching the dichroic beam splitter. The optical arrangement lor sys- 
tem II was identical to system I after the dichroic beam splitter. Additionally, a quad- 
rant star tracker with an f/3 lens was located in the elevation axis. 

Figure 3-5 shows the exterior of the elevation axis. 

The movable mirror varied its position from 4 cm to 100 cm with respect to the fixed 
system. The movable mirror traveled on precision stainless steel rods using 
Thomson ball bushings to reduce frictional loading, and was driven by a gear head 
servo motor. The location of the movable mirror was controlled by an electrically 
operated system with readouts at tho ground station to indicate mirror positions. 

The electronic system for the optical receivers was designed to extract the atmos- 
pheric modulation on the laser beam and convert it to an analog signal of 1-UHz band- 
width. Since It was desirable to observe low-level signals accurately, it was neces- 
sary to eliminate the effects of background. Therefore, the laser beam was chopped 
with a 50% on-off duty cycle at the transmitter at a 5-kHz rate and transmitted through 
the atmosphere to the receiver, thus acquiring the atmosphere modulation. The h-.ser 
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Figure 3-4. Optical Schematic, Elevation Axis and Optical Detectors 
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receiver observed both the chopped laser beam (with atmospheric modulation) and a 
slowly varying background and rejected the background with a high-pass filler. Th:. 
signal was returned to the ground (via the telemetry system). Figure 3-G ; hows the 
oplicul receiver electronic system. 

li. 1.3 QL'ADHAi, T STAR THACKElt 

A quadrant star tracker mounted on the elevation axis was used to track the 0.51-15- pm 
laser beam from the ground transmitter to an accuracy of *1 milliradlan. This infor- 
mation was relayed to the altitude stabilization system to point the elevation axis with 
its four optical receivers. Since the lens in front of the quadrant photodcleelor made 
such a small image that the range of the star tracker was leas than that required, it 
was defocussed to provide n greater range of error signals. 

The quadrant axes of the photodetector were positioned -15 degrees (see Figure 3-7) 
from the elevation and azimuth tuxes so that the error signal for each axis generated by 
differentially connecting the diagonal quadrants was 

AA 2 J 2 r Ax for Ax « r 
for uniform energy density, where 

A A ■* increase in spot area in quadrant 1 minus decrease in spot 
area in quadrant 3 

Ax » misalignment along x axis 

r * spot radius 

Positioning the Instrument axes -15 degrees from the quadrant detector axes resulted 
In no crosstalk between the error signals for small misalignments. 

The 1-lnch useful pholocnlhode surface allowed acquisition of the ground transmitter 
beam over a s 12. 5-degrec field of view. 
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3.1.4 ATTITUDE STABILIZATION SYSTEM 


To keep the optical receivers continuously pointing at the upcoming beam, it was 
necessary to mount them within two controlled orthogonal axes of an elevation-under- 
nzlmuth configuration. These axes were controlled fiom the ground by way of the up 
telemetry link, and their positions were known at all titr.es from data transmitted on 
the down telemetry link and displayed for the operator on the ground. Once the up- 
coming beam was acquired by the quadrant star tracker, it developed azimuth and 
elevation error signals which were amplified and applied to the axial drives in such a 
way as to maintain the I islrumem pointing at the beam. The positioning system is 
shown in Figure 3-8. 

To acquire the beam initially, the opernlor utilized the elevation and true bearing data 
from the beam pointing equipment. lie first commanded the elevation axis to depress 
by the same angle that the 0. 5145-pin laser beam from the ground tracking system was 
elevated. Next, lie commanded azin ith axis of the payload to look along a true bearing 
that was the reciprocal of the bearing front the ground tracking system to the balloon. 
When beam interception did not occur, the oper- tor conducted a search through a small 
zone about litis nominal position In effect intercept. An auto gale circuit automatically 
put the attitude stabilization system into the automaiic trncking mode when the beam 
was detected. An override command was provided in order to break track on any un- 
dcsired energy sources. 

The relationship of the gondola to magnetic north v.as measured with a flux gate trans- 
mitter of the type used in gating aircraft directional gyroscopes. For best accuracy 
and consistent operation, tho device was located as far- away from the rest of the equip- 
ment as possible. For this reason, it was placed at an outboard location on one of the 
gondola arms. Thus, the flux gate transmitter measured the relationship of the 
gondola with respect to the earth. The optical receiver syolem was positioned in azi- 
muth with respect to the gondola and the angle between the two was measured by a syn- 
chro differential transformer which was driven electrically from the flux gate trans- 
mitter. The output of the synchro differential was then the bearing of the optical re- 
ceiver system with respect to magnetic north. This output was converted to a dc 
analog voltage which fed a voltage-controlled oscillator (VCO) in the telemetry system. 
A servo re(»eater shown in Figure 3-8 accomplished this function. 



Each axis was driven by a dc servo motor and amplifier combination. The amplifier 
inputs wete cither the amplified signals from the quadrant star tracker in the automatic 
tracking mode or a constant velocity (dc level) input that was applied when manual con- 
trol was lieing exerted from the ground station. The groi. d controller could also 
override all other controls by actuating the standby mode. This mode was employed 
to break track from any undesired energy sources. At any time that control was not 
being exerted on the axes, they wete held in position by eieetro-mecnanieal energize- 
to-release type brakes. 

When the optical receiver system was positioned bv either axis, the remainder of the 
equipment (includ'ng the balloon) reacted in the opposite direction. The amount of the 
reaction was proportional to the ratio of the mom- nts of inertia of the opt al receiver 
system and the remainder of the equipment. When the equipment reacted to an eleva- 
tion change, the rc-olur lug fcrc- of gravity quickly restored the original orientation. 

The restoring feces in azimuth were less dependable -aid of lesser magnitude, how- 
ever, consisting of wind effects and friction between the equipment and the surrounding 
air. To keen the apparatus relatively stable then, it was important to minimize the 
reaction to azimuth motion. This was accomplished by utilizing the gondola ns a re- 
action wheel. As little equipment as possible and practical was contained within the 
positionable portion of the payload, and the remainder was packaged al the extremities 
of the gondola. In this manner it was possible to keep the ratio of the azimuthal mo- 
ment of inertia of the parte below the azimuth axis lower than one-tenth that of the re- 
mainder of the equipment, a condition llu.t assured relatively smooth operation. Thc- 
arms of the gondola also served the purpose of protecting the centrally loomed optical 
receiver system (the most expensive part of the equipment) upon landing. 

3.1.5 GONDOLA 

Tho gondola (Figure 3-9) was the basic framework for supporting the remainder of the 
systems of the balloon payload. It consisted of an aluminum girder frame supporting 
tho azimuth axis and slip ring at its center, and the three battery boxes at its extremi- 
ties. In addition to its role a3 a framework and a container for the batteries, the 
gondola also acted as a reaction wheel for the servo system and as a protective cage 
for the more delicate systems of the payload. Because of its function as a protective 
cage, it normally sustained some damage in each flight and was therefore constructed 
as simply as possible to facilitate repair. Corrcgnled cardboard crush pads below 
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each battery box absorbed most of the landing Impact and were replaced after each 
flight. The top of the gondola was equipped with a roll cage to protect the azimuth 
servo in the event of a rollover, and with sheet metal covers to act os sun shades as 
well as deflectors for ballast dropped from the balloon. 

The gondola was basically a hexagonal paralleloplpjd with a height of approximately 
5 loot and width across the flats of the hexagon of slightly less than 10 feet. It weighed 
(Ineluding the parts of the azimuth and elevation axes) approximately V50 lbs. Bat- 
teries were capable of operating the payload for 6 hours. The gondola was suspended 
in flight by nine steel cables (three from each arm) which wore fastened to a steel 
O-ring approximately 15 feet above the azimuth nxis (see Figure 3-1). 

"'his method of s- pport eliminates the bending loads which would be incurred by a two- 
or three-point support, and prolocts the payload from damage by the shock of the para- 
chute open! ig. 

3.2 THKBMO.SONDE 

The optical measurements performed have value only for the ad hoc conditions present 
during the experiments unless the atmospheric conditions are known well enough to 
allow a fu'l understanding of the turbulence effects and provide a real comparison with 
vertical propagation theory. This implies that not only the macro-meteorological 
variables such as wind speed and temperature as a function of altitude must be known, 
but that micro-meteorological variables such as the fluctuations in the temperature 
must also be known. The purpose of this section is to discuss the measurement basis, 
techniques, and apparatus for the fast temperature fluctuations. 

For "clear” weather, the single parameter that describes the atmosphere's adverse 

effects on light propagation Is the structure parameter of the index of rofr .ction*. 

This structure parameter characterizes the Index of refraction fluctuations and enters 

1-5 

into the theoretical description of the atmosphere's effects o.i light propagation ; that 
Is, variance and covariance of scintillations, aperture averaging, transmitter aperture 
dependence, beam wander, beam spread, and the phase-iront effects. 

Being a property of a gas (air), the index of refraction fluctuations can depend ou Its 
constituents, 11s pressure, ana its temperature Experi mental ly , it has boen shown 

a* 

that the only dependence on the constltutuents is via the density , which is well known. 
In addition, the pressure fluctuations sro negligible 0 relative to those of the tempera- 
ture. Therefore, the temperature fluctuations are directly related to those of the tndox 
of refraction. 


3-16 





I 


n 

Formally, the Obukhov-Kolmogorov theory of turbulence states that the structure 
function of the index of refraction behaves as 


iy P ,h) = < (n (Xj) - N(x 2 ))S-. c N 2 (h) p 2/3 <D 

2 ' 
for P = x^ - x 2 and the index of refraction structui o parameter, where h is the 

altitude, and is valid only with l o «p«L 0 , and liio line joining the measurement points 
is perpendicular to the wind direction. 


In like manner, the temperature structure function 7 j 8 
and 7 


»>•!-(/>. h) C. r 2 (h) p 2/3 




(2> 

(3> 


for pressure (P) in millibars, temperature (T) In degrees K. , and wavclei gth (\) in pm. 
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A measurement technique that suggests itself (and is widely used in ground systems ) 

is to mount two fast, sensitive sensors in a bridge configuration at h known distance 

opart. My suitably squaring and averaging the output, the temperature structure func- 

•> 

tion is obtained. With p known, equation (2) is used to obloin C_ ,** (h) and equation (3) 
.2 t> 

is used to obtain (h), using density tables and a simultaneous measurement of the 
temperature. 

Sensitivity, speed, and uniformity of tho two fast sensors was of prime importn eo. 

O 

High signal levels were required to accurately measure D T " and the sensors had to 

have matched temperature-resistance curves over the entire range of temperatures to 
14 

be measured . With thin-wlre sensors, speed was no problem, and uniformity of 
length and cross-section was relatively easy to obtain. Sensitivity was the main ques- 
tionable feature of ihin-wire sensors. 

"The responslvity of a metallic wire sensor can be calculated from the applicable re- 
sistanec-vorsus -temperature equation 

Rj. “ R 0 (1 ♦ a T) (4) 

where R q = p//s(pboing the roslstlvity, l tho longth, and s the cross-section) and a is 
the tomperature-resistnnee coefficient. 


V 
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The value p for plat ium-rhodium^' is of order 20 x 10 1(> ohms per cm, and a is of 

i o __o q 

order' *.3x10 per C. Therefore, lor a 1-nwn long. 2- fi diameter plntinum- 

rhodium wire, R ~ 0411 andAHft 8,3 x 10 “ ohms per °C. Therefore, domi.r.ding a 
o , • 

sensitivity of 5 x 10 °C means a resistance difference of 41.5 pohms must be dis- 
cernible. 

Too r.u.eh current cannot be fed through the wire Iwcause the self-healing will mask 
and distort, the tuitadence effects the wire *s sup|>oscd to measure. 

Therefore, with these design restrictions nd the fact thnt a temperature difference 
was desired and amplifiers that do not extet.\ down to dc have much lower noise, the 
circuit used is as shown* 1 * in the block diagram. Figure 3-10. 

The thermosonde is the box shown or. the right hand sidi of Figure 3-13 susjx .ided 
from the load bar. The rather large buik Is due to .he large bailer) supply recji red 
to run the system for 8 hours and to keep It warm for that period. The actual circuitry 
a. id instrumc Ta'le". weighed less than 20 lbs. In order to eliminate the effects of 
balloon wake in the C^, measurement, the llicrmoscnde was -uspended from a reel 
hung from the load bar of the balloon. When the balloon reached an altitude of 330 m 
(1000 ft), the reel was released allowing the thermosonde to descend. A hydraulic 
drag mechanism in the reel limited the rate of descent io IV ms * (30 ft s *) giving n 
deployment time of approximately 20 sec. This deployment technique allowed error- 
free measurements to start at an altitude of about 60 m (ss200 ft). The reel mechanism 
and the sonde can lx? seen in Figure 3-14, 

3.3 FLIGHT VEHICLE 

3.3.1 BALLOON CHARACTERISTICS 

The balloon selected for the BAPE flights was a 2.0 x 10 ( ’ ft 3 polyethelene balloon. 

This selection was based upon a requirement of an altitude of about 90,000 feet (27.4 
km) and the gross weight of the BAPE payload (»750 lbs) plus the weight of Air Force 
balloon control and navigation equipment (ftslOOO lbs). The total weight of approximate- 
ly 1750 lbs plus the iltitude of 90,000 feet makes the 2 10*’ ft^ balloon the proper 

choice for this type of flight (see Figure 3-11). Launch, navigation, recovery and 
all other assets of tie flight vehicle were handled by Detachment 1, Air Force 
Cambridge Research Laboratories — Balloon Research and Development Test Branch , 
Holloman Air Force Base, New Mexico, under the command of 1.1 Col Hobort J. 

Reddy. 
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Figure 3-11. Balloon Requircmen.j 



3.3.2 LOAD BAR AND RIGGING 


The expel iment flight configuration ia shown in Figure 3-12. The BAPE payload Is 
supported by nine steel cables which converge at a central point (Figure 3-13). This 
l»oint is attached by an explosive shackle lo the Air Force load bar which supports the 
Air Force navigation and flight control system as well as the thormosonde system and 
reel, in Figure 3-13 the components supported on the load bar arc, from left to right: 
ballast hopper, command receiver and telemetry system, backup navigation transmit- 
ter, ba’last hopper, and reel and thermosonde assembly. Flight configuration (reel 
not released) is shown in Figures 3-14 and 3-13*. The load bar, as shown in these 
figures, Is supported by a second set of steel cables which converge on a. ring at the 
bottom ot the parachute shrouds. The ioad is supported on the shrouds of the open 
parachute and attached to the bottom of the balloon. The balloon is of the zero- 
pressure type meaning thal the bottom of the oalloon is oocn and the helium ia trapped 
only by its own buoyancy. 

Vertical control of the balloon is controlled by the ballast l;opf>cr, shown in Figure 
3-13, which drops iron filings to ascend, and n poppet valve nooui Its inches in diame- 
ter (Figure 3- id) which releases helium to descend. Horizontal control of the balloon 
is obtained by ascending or descending into wind layers blowing la the desired direc- 
tion. 

I’poa termination of the flight, ihc balloon is destroyed by nn explosive device and 'he 
puyload descends on the parachute (Figure 3-17). Upon impact ttie chute and load bar 
arc disengaged from the payload by the explosive shackle so 'hat ground winds do not 
drag Hie payload. 

3.3.3 EXPERIMENT KITE AND FLIGHT PROFILE 

The experiment was performed at the SCAT site on Holloman Air Force Base which is 
pnrt of the White Sands Missile Range In New Mexico. Tho position of-lhe silo is shown 
tn Figure 3-18. The site was on a level plain with semldeaurt vegetation approximate- 
ly 7 miles northwest of Alamogordo, New Mexico. The plain is a valley floor with an 
altitude of 4000 ft (1.2 km) above sea level. The San Andres mountain range with 
peaks of about 8000 feet (2.1 km) lies approximately 30 miles to the west, and the 
Sacramento range with peaks to 11,000 fcoi (3.3 km) lies 10 miles to the east. During 
the early morning hours, cold air trapped in tho valley during the night drains south- 
ward. 

*T‘<c balloon in Figure 3-15 is a 60— ft -diameter balloon much smaller than tho balloon 
used for the flights. It was used for tethered teste. 
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Figure 3-12* Flight CoafiKurniion 
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Figure 3-15. HAl’F l-'light CuiiiigiuiUiuii 






Figure 3-16. Poppet Valve 
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Tho balloon was launched at 0G13Z (12:15 AM MST) from a point near Upham, New 
Mexico (50 miles went of the tracking site) and slowly drifted eastward arriving over 
the tracking site at about 10-1 tZ (4:-H AM MST). Tracking commenced at about 113CZ 
(5:.')0 AM MST) and was completed .it 1415Z (8:15 AM MST). The horizontal flight pro- 
file is shown in Figures 3-19 and 3-20. Altitude, . enith angle, and slant range are 
shown In Figures 3-21 and 3 -22 and tabulated in Table 3-1. 

3 • •* GROUND TRACKING STATION 

The following sections dcscrilo the mobile laser tracking station which was used to 
(rack the balloon-borne payload, with laser beams. The system cons'sted of three 
basic units which were transported to the : te in three trailers w ee Figure 3-23). 

An artists conception of the system Is shown in Figure 3-24. The right-hand 
trailer as shown in the cutaway contains the lasers, modulators power supplies and 
other equipment necessary to operate the laser beams. The output is transmitted 
through a series of mirrors in the tracking mou towards the target. The second 
trailer contains the tape recorders, balloon control system telemetry, and other in- 
strumentation necessa. y for the experiment. In actual operation the main optics of 
the lei- scope were used only for reception purposes, and the laser beams were trans- 
mitted from a |>ort to the right i f the main telescope. The artists conception is in- 
correct In this aspect. The system was developed for the BATH experiment at GSFC 
Code 524. 

3.4.1 LASER TRANSMITTER SYSTEM 

The laser transmitter system, which was housed in tne instrumentation trailer shown 
in Figure 3-23, projected its output through a tube into the tracking mount where a 
scries of mirrors directed the radiation to the target. The transmitter system, as 
shown in Figure 3-25, consisted of two lasors: an RCA Model J-15268 and a Sylvania 
Model 948-1. Table 3-2 lists the characteristics for each of the lasers. 

The lasers, wldch were bolted to ^ steel table Inside the instrumentation van, re- 
mained stationary throughout the experiment. The laser beams , after being passed 
through collimation and alignment optics, were projected coaxially into the tracking 
mount. The output characteristics for the transmitter are shown in Table 3-3. 
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Figure 3-19. Vertical Flight Profile 
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Tnble 3-1 

Pnrnmvterg of Bnlloon Trajectory, Uciol>cr 21, 1970 
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Figure 3-24. BAPE Mobile Tracking Station 






















Table 3-2 

Laser Characteristics 


Characteristic 

Value 

RCA Model J-15208 

Sylvania Model 948-1 

Type 

Argon ion 

CO2 molecular 

Mode 

CW 

CW 

Wavelength (microns) 

0.5145 

10.0 

Maximum output (watts) 

4.3 

15 

Bourn divergence (milliradians) 

2.0 

5 

Beam diameter (millimeters) 

4 

6 


Table 3-3 

Transmitter Characteristics 


Characteristic 

Value 

Argon 

CO2 




Wavelength (microns) 

0.5145 

10.0 

Mode (sine wave modulated) 

5 kHz 

5 kHz 

Power out (mean watts) 

1.0 

5 

(peak watts) 

2. 0 

10 

Divergence (milliradians) 

1.0 (1/e") 

1.0 (1/e 2 ) 

Ncnr-ficld diameter (millimeters) 

0 (l/e 2 , 

30 (1 /«»2) 


The COy and argon laser beams, as shown in Figure 3-25, passed through a 
diehroic beam splitter, an attenuator, a 5-kIlz modulator, and a collimator (5x for ‘he 
CO2 beam, 2x for the argon beam) to diehroic beam combining mirror. The com- 
bined beams then passed through a safely shutter to the two mirrors that aligned the 
beams to the tracking mount. The first of these mirrors was slightly transparent in 
the visible range, which allowed a target alignment scope to view a direction parallel 
and coincident to the outgoing laser beams. A beam direction monitor, located between 
the alignment mirrors, allowed alignment of the target allgament scope to the laser 
beams, but was removed when the system was in operation. The targe - alignment 
scope Imd a protective filter for the eyepiece which allowed the operator to view (he 
target during laser operation. During tracking, the operator could view the target in 
the alignment scope through the coelostat of the tracking mount to ensure that the laser 
beams were on target. The operator adjusted the alignment scope crosshairs onto the 
target by means of the last alignment mirror in the transmitter system. 

Power monitors, mounted in the system, monitored the output power of the lasers and 
the modulation frequency. Because of the difficulty and danger involved in aligning the 



transmitter optical components with COg and argon lasers, a helium-neon laser, also 
mounted on the transmitter table, was used as an alignment tool for the optics. 

3.4.2 TRACKING MOUNT 

The tracking mount, shown in Figure 3-26, received the radiation from the laser 
transmitter system and directed the radiation through a coelostat system to the target. 
A diagram of the coelostat system is shown in Figure 3-27. Radiation from the lasers 
traveled horizontally through a pipe from the laser van to a mirror at the bottom of the 
azimuth axis.. (The pipe from the van to the tracking mount is shown in Figure 3-24.) 
This mirror deflected the laser beams to a second mirror at the top of the azimuth 
axis, which then defected the beams to a third mirror at the end of the elevation axis. 
The third mirror deflected the laser beams to a fourth mirror located slightly to the 
right of the main telescope. This mirror directed the laser beams through a small 
lube which is, as shown in Figure 3-20, to the right and parallel to the line of sight of 
the main telescope. The laser beams were directed through this tube Without being 
collimated by the receiver optics. 

The receiver optics, a 450-inch (11.4 m) focal length, 30-inch (70 cm) diameter ens- 
seprain, were used to collect laser rudiatio . reflected from the cube corners mounted 
in the balloon payload, and to f oous the radiation on a star tracker at the back c.f the 
telescope (see Figure 3-20). The collected radiation was used by the star tracker as 
a signal for automatic tracking of the target. 

The tracking accuracy of the star tracker was approximately 5 arc seconds rms; 
however, because of intermittent loss of autotrack due to severe scintillation during 
the experiment, the target had to be manually tracked by the observer (shown in Fig- 
ure 3-28). Manual tracking was accomplished by observing the target through' a guide 
telescope and adjusting .he main telescope by means of a Joy-stick control mounted 
on the manual control panel. Details for the guide scope and control panel arc shown 
In Figures 0-26 and 3-28. Tracking accuracy in the manual mode was approximately 
30 arc seconds rms. 

In the autotrack mode, the star tracker intermittently lost track and went Into the 
search mode because of insufficient integration time In the alar tracker circuitry. 

After this experiment, modifications were made to correct ihls problem with the star 
tracker circuitry, and the star tracker tracked perfectly for BAPE II flights during 
September of 1071. 
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Figure 3-26. BAPE Tracking Mount 
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The tracking mount is a converted NIKE -AJAX tracking pedestal to which the 30 inch 
telescope and eoelostat system have been added. During operation, the tracking mount 
is elevated by the aluminum legs shown in Figure 3-24 so that the elevation axis is 
15 feet (5 meters) aSx>ve ground level. The tracking mount is elevated for the following 
reasons: 

To give a clear view of the horizon over trnilers ard buildings in the 
surrounding area 

• To provide a safety factor by teeping the laser beams well above 
personnel in the nrca 

• To keep the laser beams above the low level ground turbulence which 
does not follow the norr nl Kolmogorov spec!' urn of turbulence 

3.4.3 TELEMETRY, PAYLOAD COMMAND, AND DATA HECOxtDIN J SYSTEM 

Command of the balloon payload and data from the balloon payload were obtained 
through an L- and S-bn<ul telemetry system. /' console (Figure 3-2!') In the instrumen- 
tation van allowed an operator on the ground to control the azimuth elevation, and 
operational stilus of the payload. Two small parabolic telemetry antennas (Figure 
3-30) transmitted the information to and received information from the payload. In 
addition to the command and telemetry system for the payload, a standard radiosonde 
transmitter operating between 1000 and 1700 Mil/, was used io ti ansinit data from the 
thermosonde to the ground. Telemetry characteristics are listed in ruble 3-1. 


Table 3—1 

Experiment Pavlcnd Command and Telemetry 


Characteristic 

Valuo 

1431-MIlz command system 

Conic Corp. Model CTM-lJHF-305 

Output 

0 watts 

Antenna 

4 -foot paraoblic, 10-inch primary 
bcnnnvidth 

2241.5-MHz telemetry 

Conic Cory. Model C FM-UHF-305 

Output 

5 watts 

Antenna 

Hemispherical pattern, omnidirectional 

IGG0- to 1700-Mllz telemetry 

Viz Mfg. Model 1001, military Model 


No. AN/AMT-12A 

Output 

400 milliwatts 

Bandwidth 

500 Hz 

Modulation 

100','c, audio 


- ;• • 
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Downlink telemetry carrier assignments arc shown in Table 3-5. The outputs of these 
channels ..ore recorded on an I'M tape recorder. Data recording assignments and 
recording circuits are shown in Table 3-6 and Figure 3-31. respectively. 


table 3-5 

Downlink Telemetry Carrier Assignments 


Channel 

Function 

3B 

' "" 

Spare 

on 

Spare 

ini 

Argon movablo 

1311 

. CO^ fixed 

711 

Spare 

911 

Argon fixed 

1511 

CO., movable 

•i 

Tracker acquisition Indicator 


Pressure 

0 

Balloon receiver ACC 

7 

Commutated 'emperntures 

8 

Mirror position 

9 

Azimuth 

10 

Elevation 

11 

Spare 

12 

Azimuth error 

13 

Elevation error 




Proportional bandwidth 



";7 
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Table C-6 

Data Recording Channels 


T.R. Channel 

Signal 

Bandwidth 

Record Module 

#1 

Argon A fixed 

1 kHz 

FM 


#2 

Argon 3 movable 

• kHz 

FM 


* 3 

CO^ A fixed 

1 kHz 

FM 


#4 

CC >2 B movable 

1 kHz 

FM 

1 

#5 

36-bit BCD time (1-kHz carrier) 


Direct 

#6 

Balloon pack az error 

DC 

FM 

. i 

#7 

Balloon pack el error 

DC 

FM 

i 

iS 
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DC 

FM 


#9 
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DC 

FM 


#10 

Ground rec argon demod 
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FM 

' 
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FM 

i 

? 
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Ground comm data 

DC 

FM- 
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#13 

Trans power out (argon) 

kHz sin) 
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DC 

FM 
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BAPE commutated temp 

DC 

FM 
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#2 
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(1-kHz carrier) 

Direct 

10 6. Balloon pack. el 1 

#3 

Grd thermal sensor sig 

; kHz. 

FM 

M 7. N/. A position az 

#4 

Crd thermal sensor sig 

1 kHz 

FM 

M 8. N/A position el 

#5 

Wind velocity 

1 kHz, 

FM 

9. Calib 

#6 

rec argon 

5 kHz ± 1 kHz 

Direct 

10. 

#7 

Grd rec CO„ 

5 kHz i 1 kHz 

Direct 




Figur'' 3-31. Data Recording Assignments and Circuits 
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4. CHRONOLOGY OF EXPERIMENT 


This section explains the time sequence of events which occurred during the experi- 
ments and gives a chronological description of the events. 

The balloon was launched at 12:13 AM MST from Upharn, New Mexico, a point nearly 
50 miles west of the tracking site. The balloon ascent rate was set at approximately 
1000 feet per minute. At approximately 12:15 AM MST whei *he balloon had risen to 
an altitude of 1000 feet, the reel supporting the thermosondc was released and about 20 
seconds later the sonde was fully deployed 1000 feet below the main payload. At 12:30 
AM MST the flashing strebe lights of tt'e balloon were spotted from the tracking station 
and from that time on until termination ■' the experiment the tracking station was in 
visunl contact with the balloon. The weather was clear with no clouds in sight during 
the experiment. High winds in the 40,000 to 80,000 ft altitude range droi • the balloon 
quickly eastward at speeds of up to 90 knots, nH by 1:00 AM Me f the payload was close 
enough to command b> tc Ir metry. A check of payload operation was made and every- 
thing was fourd to be in order. During ascent, a radiosonde receiver at the launch 
site was able to receive data from the thermosonde to an altitude of about 20,000 feet. 

A soennd receiver, : ocated at (lie tracking site with a higher gnin antenna was sup- 
posed to pick up the data collection at this point but due to a jnalfunctlon was unable 
to operate. A b° roswiteh on the thermosondc automatically shut it off at 45,000 feet. 
After reaching a float altitude of about 92.000 iect, the balloon was slowly manuevered 
towards the tracking station by slightly changing altitude to allow the prevailing winds 
to drive the balloon towards the tracking station. Because oi the low velocity of the 
winds, the Ir.Uoon was not within tracking range until about 5:30 AM. At 5:37 AM 
tracking and data taking was started and continued until 6:10 AM, at which time, data 
taking was stopped in order to make a quick analysis of the data to make sure that all 
systems wore working correctly. At 7:05 AM tracking was resumed and continued 
until the balloon went out of range at 8:00 AM. Dawn occurred at 7:13 AM, nnd there- 
fore most of the data was taken during the dawn calm, the period of lowest turbulence. 

At about 7:00 AM a 3low descent of th-- balloon was started as shown in paragraph 3.3.3 
and at 8:30 AM the flight was terminated. The payload descended by parachute from an 
altitude of 45,000 feet and landed near Rosw«. U, New Mexico. 
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5. DATA REDUCTION 

Tho following sections discuss the method of data reduction used in reducing 'he experi- 
ment data from raw form to final results. 

5.1 THERMOSONDE DATA 

Table 5-1 Hr ■ i U. > strength of turbulence versus height as obtained from the therm i- 
sonde. Figure 5-1 allows the product of strength of turbulence, and height to the 5/6 
power versus height. This product controls the strength of scintillation and therefore 
tho figure shows that ihc bulk of the scintillation caused by turbulence below 6 km was 
caused by a layer at approximately 700 meters above the ground. 



Figure 5-1 
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Table 5-1 
Thermosonde Data 


h 

| (h-hR) 

<Mt R ) 5/C 

(h-hg) 5//3 

c 4 

C N 2 (h) (h-h R )° /6 

0 N 2 <h) (h-hg) 5/3 





L’nfts 



m 

m 

,„ 5/6 

5/3 

m 

-2/3 

m 

1/0 

m 

" 

m 

1-190 

290 

113 

1 . 29 x 10 * 

-14 

0.18 7. 10 

20.3 x 10 14 

2.30 X 10' U 

1560 

300 

134 

1.80-8 10 1 

0.19 X 10“ 14 

25.0 X 10“ 14 

3.4. X 10 -11 

1 (550 

-150 

101 

2.09x10“* 

t. 19 X 10" 14 

31.2 X JO* 14 

5.11 X 1J' U 

1850 

(150 

2i 0 

4.84 X 10 4 

0. 04 X lo' 14 

141 X 10 -1 ' 1 

3. 10 x 10~ 10 

19 SO 

"80 

260 

0. 70 X 10 4 

0.20 X 10' 14 

(57.8 x lo“ 14 

1.70 > 10 -1 ° 

2 ~0 

1270 

390 

1.44x10° 

0.39 x lo" 1 *’ 

1.48 x JO" 14 

5.03 20“' 2 

2590 

1390 

41) 

1 . G8 X 10° 

0.24 ;< 10' 1G 

0.98 x 10 _1 4 

4.03 x 10' 12 

2900 

1700 

500 

2.50X10° 

0. lb \ 10' i0 

0.90 x lO -1 * 

4. 50 J. 10 _i * 

3-100 

2200 

500 

3. 00 x 10° 

0.12 x 10' 1(i 

0.72 x 10* 14 

-12 

4.32 x 10 

37-10 

2510 

6 Hi 

4.02 x 10° 

0.15 x 10 _1 ° 

1.02 x 10’ 14 

6.9* X 10~ 12 

5070 

38/0 

9 SO 

9. 00 X 10° 

0.21 x 10' 5 ° 

2.00 x 10' 14 

2. 92 x IP -11 

5 115 

3915 

1000 

1 . 00 x 1 0° 

<>. 17 x 10~ 1G 

1.70 x lo’ 14 

5.70 X 10' U 

5185 

•1295 

1050 

1. 12 X 10*’ 

0 . 10 x 30 10 

1.00 x lo' 14 

1.12 X 10" U 

5070 

4-170 

1100 

1.21 X 10 U 

1.2 x 10' 1G 

13.2 x 10' 14 

1.45 X 10* 10 

uH7i) 

4(575 

1130 

1.28 X 10° 

0.079X lO' 1 *’ 

0.80 X 10" 11 

1.01 x !0 -J1 

0950 

4750 

1140 

1.30x10° 

— IP 

0.41 X 10 

4.0(3 X 10 14 

,,-11 
5.33 x 10 

0090 

4890 

1190 

1.42x10° 

0. 097 x 10' 1G 

1.15 x 10~ 14 

1.37 x 10 _U 
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5.2 


METHOD OF OPTICAL DATA REDUCTION 


Data reduction for the two channels of argon detector data and iho two channels of C0 0 
detector data yielded the important statistics of probability densitv function, cumulative 
density function, variance, autocorrelation, power spectral density, and cross-cor- 
relation. With the exception of power sp-ciral density, all statistics were obtained by 
playing back the tape-recorded analog data into a Hewlett .’nekard Model 3721A cor- 
relator. This instrument is a digital signal analyzer which derives Ihe desired stntir- 
tics by sampling the input data and operating on the digital result with hardwired sta- 
tistical functions. All statistics wc.-e calculated ir the summation mode in which N 
regularly spaced samples of data arc used to compute the desired statistic and the 
result is divided by N. The sampling rate was set at 10'* sec * for the argon detector 
data and 10" sec 1 for the CO,^ detector data. IT. sc were chosen froi • estimates of 
signal bandwidths to be small enough so that io aliasing crro.s were pres' nl. The 
value of N was set at 3.2 x 10* or G.4 x 10 which provides, for example, a statistical 
average over 32 or b-4 seconds of data when a 0. 001-second sampling interval is used. 
The resultant statistic was displayed as an X-Y plot on i CItT display or hnrd-copied 
on x standard X-Y ink recorder. The output display hud a resolution of loo points 
equally spaced across the horizontal axis. 

Resolution wi.s 200 levels lor the full vertical axis nnd was complemented by inter- 
polation between points on the display. 

The first processing of the data was done in the autocorrelation mode. First, the dc 
level of tiie input analog data was removed and then the autocorrelation was computed. 
Argon detector data analysis employed those coupling capability of the correlator with Us 
1-Hz cutoff frequency. CO., detector data was first passed through an external high- 
pass filter with a 0. 1 cutoff frequency nnd then Into the dc -coup led correlator. This 
was done to preserve more of the typically low frequency CO,, detector data. 

The value of the autocorrelation function at time Ing zero (mean square value! was 
read from the CRT display and written down br a function of detector number. Simul- 
taneous data segments were analyzed for each pair of detectors. Then the cross-cor- 
relation between like detectors was computed and its value at time lag zero recorded. 
This value divided by the square root of the product of both individual mean square 
valuen from the autocorrelation analysis was the desired correlation coefficient. All 
scale factors between input irrndloncc at the detector and tape recorder output voltage 
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cancelled out when the coefficient was formed. The correlation coefficient was com- 
puted for several detector separations and the correlation function for each wave- 
length was determined. 

The mean square values obtained above were just the data variances since the dc or 
average value was zero in each case. The variances were in terms of tape recorder 
output voltage squared. It was desirable to normalize the variances by the appropriate 
dc values squared in order to have a standard summary statistic for the strength of 
irradiance fluctuations. The dc value for each data record was obtained irom the 
SO'A level of the cumulative density function ns the data was once again played back 
into the <orrelator (this time, dc coupled). It was also necessary to compute the dc 
value of the background (signal from detector in absence of laser irrad'ance) and sub- 
tract this from the data dc vnluc before the final ratio was computed. Several cumu- 
lative dt nsity functions were nls^. plotted on the X-Y recorder and then transferred to 
probability paper to test for log-normality. 

The autocorrelation function was transferred to the X-Y recorder in the same manner 
as above and the 100 function points were fed into n digital computer. The computer 
was programmed to perform a Fourier transformation of each input autocorrelation 
function. The result was the desired power spectral density function. 

DefiniUjns of die terms used in reducing the optical data are listed below: 

N — Sample time in seconds 

p — Separation of sensers - cm 

Rj(0) - (Sj - Sj)2 

r 2 (0) -<s 2 -s 2 > 2 































OPTICAL DATA POINTS (Cont) 
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0. DATA ANALYSIS 


0 . 1 V K NT IC AL IWIIT OF Tl’IUU'I.KNCK 

Most of the estimates of the vertical profile of turlmlcnce have lieett based upon mea- 
surements of stellar scintilla (Ion and a.-di nnomlcal seeing. Hufnnjtel and Stan lev' 

•» 

proposed one of the first and be^t-known profile* and in It) ft ft Hufnajtel" revised the 
profile (Figure -1). Fried 3 projiojted n mathematical model ef the profile of atmos- 
pheric turbulence and has uses! It extensively In his numerous papers on optical propa- 
gation (Figure 6-21. 

Due to the presence of a sharp temperature gradient at the tropo pause, most authors, 
have assumed that the tro|>o|>ause Is the major source ol stellar scintillation effects. 
Analysis by Tatars!;! based upon measurer. .ents of aperture averaging bv I’rotheroe* 
and nieFtr covariance measurements by Keller 5 tern! to .oibsiantiate this hvjs'thcsls. 
Tltu'rto -' and HtooMtcr' have politic* I out that Fried'# rnfclel of the atmosphere over- 
estimates the strength of turbulence In the firs! few kilometers of the atmosphere wl,;tn 
nol properly nreoiintlny foe .‘(irt.ilrmn near (lye I rri.w\n;i,ie.i . The major it>.(nts of evidence 

Indicating a n- rrow band of high htrliulenre at or near the tn'pofiause are the sU*> of the 
correlation length ol stellar scintillation and ihe eorrel itlon of (he po *er «|K*ctral 
density of scintillation to wind sjieixi at the trnpnpausc. However, wide variations In 
the «Ue of the stellaV scintillation correlation distance (from it (o 10 cm) make the 
determination ol the altitude of the scintillation . .eneratmg latent of tnrUdenr. dif- 
ficult . 

The strength of turbulence Is normally specified by the refractive Index -st rue lure 
constant 


*(N. -N,) 2 > 


I <0' L 
o o 


where N and N., are the refractive Indices at two (mints separated bv a distance p. 

1 u 

Here wo have assumed the Kcd mogorov similarity theory of turlmlcnce to bold (or 
separation (pi between the Inner (f 0 ) and outer scale (l- rj l. The refractive Index of air 
at optical frequencies may be expressed as a function of Ihe properties of a pas and 
the wavelength through the relation 


UiUi [, . JL™15 ] , 


fl-i 






0 (Om 100m Ikm 10km I OOkm 
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ALTITUDE- METER: 


Figure 0-2. Fried’s Model 


w!u ro P Is pressure In mill fours, T is absolute temperature in °K, and X Is wavelength 
in pm. Pressure fluctuations are tnpldly damped out and, therefore, refractive Index 
changes arc mainly due to mlerotliernial fluctuations. 


Defining the temperature- structure constant ns 

2 d t ( P ) <(v T 2 r> 

C T "7^ 7273 — 


l o < 0 < L o (3) 


wo seo that it can bo easily converted to C 


N 


2 f 77. OP f. 0,007 5 . r 


,-12 


(4a) 
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liquations 3 and 4a provide the most direct method by which to obtain the vertical pro- 
file of turbulence. The temperature structure constant can be measured with a pair of 
low-incrtla- resistance thermometers which are part of a Wheatstone bridge. The out- 
put voltage of the bridge, which Is proportional to the temperature difference between 
the thermometers, may be amplified and the mean square obtained with suitable elec- 
tronic circuits. D-j,(p) combined with data on mean pressure and temperature may then 
bo used to compute C^.". Krechmer * *, lata rsk 1 1 and Lawrence*** have used this 

method to obtain data in the first few kilometers of the atmosphere employing both air- 
craft and tethered balloons to elevato their inslrumcnts. Unfortunately the data ex- 
tends only up to 3.5 km. 

f t order to compare the results of the optical measurements with theory, an Instrument 

to measure the veitical profile of turbulence was constructed and flown ns part of the 

balloon payload*. The d« vice consists of a pair of low-tnertta-rcslstance thermomet"rs 

2 

•n a Whoatsene bridge with suitable processing electronics to derive C^.“. The device 
was attached to a standard MO-31 7A/ANIT-J2A radiosonde which was suspended on n 
300-m eatitc below the main payload to avoid the turbulence effects caused by the rising 
balloon. Tlie output of the resistor thermometer system was connected to one of the data 
channels of the radiosonde and data w s tclcmetenv. to the ground In a standard telem- 
etry format. Mean pressure and temperature were obtained from the standard radio- 

'J , 

sonde sensors**. The profile of O^, obtained It shown In Figure 6-3. 

Ik'c.'uise of the rapid ascent of the balloon (1000 ft/m In or m/sec) the radiosonde pro- 
vided only a few data points within the first few Inuutredj'U'ters. Therefore, 

1 ^ 

TuUirski’s method of mean temperature* was used to supplement the thermosonde 
data. From Figure 16 of Kcfcrenco 10 wo huve approximated Tatarskl a data by the 
following equations: 

Stable Stratification (Temperature Inversion - Increase of mean temporaturo 
with height) 


C T ‘ « 4 x 10 


-2 | j.2/3 h -1/3 


T, 


( 6 ) 


Unstable Stratification (Decrease of mean temperature with height) 


C T 2 * 5. 76 K 4/3 h ~ 2/3 T, 2 


( 6 ) 


•Sec Section 3.2 


**Seo Appendix A 
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C T 2 = 2.17 x 10" 2 
C N 2 - 1.47 x 10 -14 
Unsta! 'c St rati fraction 

? _o /i o 

C t '-’ 1 . 73 * T«“ ( 10 ) 

C^, 2 <■ 1.17 x 10~ 12 z ~ 2/3 T, 2 (P « 87Gmb, T “ 280°K) (1 1) 

The experiment wrtH j»er*ort..ed und'u* conditions of stnt. o strntifrnction with an ln- 

• etslon at about iso in, as 8?*o\vn by tb» vortical da shod lino of Figure 0-3. Therefore, 

•» 

•Xi.i^tko 5 was used U> calculate C^. up to nn altitude of alxiut h0 m. 'Pie results of 
Mi'* snaly. *s «ro shown In Table 0-1 and as crosses in Figure 0-3* 


Table G-l 



! log h 

h -l/3 

T.(°K) 

O O /‘I 

C T “ 0 <m" ) 

c ^ 

I If 

1.00 

0. 4 G-l 

1.40 

1.47 x 10 -2 

1.C0 x 10 ' 14 

( 20 

1.30 

0.3G8 

2.92 

2.31 x 10 “ 2 

1.59 x 10 ~ 14 

I 40 

1.00 

0.293 

5.92 

3.77 x 10 “ 2 

-14 

2.M x 10 

1 CO 

1.78 

0.255 

10.3 

5.70 x 1C -2 

3.83 x 10 ~ 14 

[ 80 

1.90 

- - 

0.231 

12.7 

6.35 x 10 -2 

4.31 x 10 “ 14 


<7) 


h- l/ 3 T. 


( 8 ) 


h‘ l/ 3 T, (P <= 87C mb, T » 280°K) ( 9 ) 


a 


Whore T. 


Tj - T. 


log 

K » 0.4 


{ i0 N_ 
- h 2 . 


Therefore 


Stable Stratifractfon 


The profile of turbulence has been approximated by a mathematical model shown by 
tho dashed line, which obeys the equation 


r 2 „ 2 . -1/3 -h/h rt 

C N “ C No h 0 ° 


( 12 ) 
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REFRACTiVE-INDEX-STRUCTURF CONSTANT AS A 
FUNCTION Of HEIGHT ABOVE GROUND LEVEL. 

1.-00 AM NET 21 OCTOBER 1970 AT HOLLOMAN AFB 
NEW MEXICO. BARE I EXPERIMENT 

> ♦ DATA POINTS OBTAINED BV TATARS XI*S 
METHOD OF MEAN TEM°£RATURES 

• DATA POINTS OBTAINED FROM THERMOSONDES 
MEASUREMENT 
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Figure 6-3 
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where 


C S(j 2 «4x 10 ' I4 m" 2/3 


h = 1000 m. 
o 


The value of 4 represents the ordinate of the model at h = li.t, although It has no 
physical significance since the refractive- Index structure constant Is not rellnb'e so 
close to the ground, 'the scale height is the break point In the monel, which for stable 
strr.t (fraction normally occurs at the temperature Inversion. 


2 

Two Integrals of C^, which will l>e of Importance In later calculations of the scintillation 
statistics are 


and 


I 

1 


C N 2 (h) dh = C No 2 h 2/3 rc/J -1.30 C Nq 2 h o 2/3 (m l/3 ) (13) 

= S.?0xl0~ 12 m l/3 

C N . 2 (h) 1. 5/G dh - C No V V2 r ( 3 / 2 ) - 0.89 c No V /2 (!•») 

t> - 1. 13 x 10~ 3 m 7/G 


Wind speed and mean temperature as n function of altitude are shown In Appendix A. 
No thermosonde data Is available above 0 kilometers; therefore, Hufnagel'r model for 
the tro|>opause turbulence la used with the integral over the tropopnusc e<|ual tc 


1. 


C 2 (h)dh - -1.31 x 10“ 13 . 1/3 

trop ^ 

Ihe final model of the atmosphere used for calculation o.’ optical propagation statistics 
l>o comes 


C 2 ■ r* ^ l" 1/3 - h /b o c rt. i v c + 

C N C No h e + 5(h-lO C No 


t 1?) 


where 


C * 
No 


4.0 x 10 


-14 -2/3 


m 


C Np " 4,3 x 10 ” 13 ml/3 


ho * 10C0 n. 


h -13.2x10 m 
P 
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fl . 2 LO O- AMPLITUDE VARIA NCE A MD SCfNTrLLATION 

The strength of scintillation Is often specified by the normal! ted power variance at V. e 
receiver; 


"Idle this is a simple parameter to measure. It Is not i. convenient term to handle 
statistically. Hecause the power Is log-normnlly distributed (t.'at Is the logarithm of 
f'ver Is normally distributed) ns a conseciuence of the central limit theorem. It Is mere 
co venlent to make calculations baaed on 



which has a not nal distribution. This ,.arni etcr, sometimes referred to ns Ior- 
power or log Intensity, may bo analyzed using tho normal distribution tables , and Its 
variance Is simnlv related to l!:e normalized power variance 

<(ln P/P) 2 > In (1 •» a ”) (2) 


Since most of the theoretic.*!! analysis of optical projmgatlon Is done In terms of phase 
and amplitude, tho analyst uses the term log-amplitude 



which also has a normal distribution. The variance of this parameter is easily related 
to the log-power (or log-lntenaUy) variance 

<(ln A/A o > 2 > ■«;* «ln P/P o ) r > -{ln(l+o p 2 ) ( 3 ) 


hi the case where 



<< 1, then 


2 

P 


(•») 


and we can make the approximation 

2 o “ 

<(ln A/A ) > « — P— 

o 4 


( 5 ) 
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Log- amplitude variance has been represented by numerous different symbols by various 
authors, but wa will use the terminology of Fried* and use C|(0) or Oj 2 . Small letter I 
will be used to represent log-amplitude. 

If n laser In-am Is transmitted from an optical system which allows for the adjustment 
of tho aperture and radius of curvature of the emergent wavefront, the scalar wave 
function at the optical system output {/. «■ 0) Is * 


O 

Ho <r)jl-' 0 « exp [~.£tT * 


where k - ?tr/X. P la the radial component of r, and t Is the axial compor.nt of r. Lot 
the beam have a central amplitude of unity and n normal radial distribution with varl- 


unce a then from Schr/oltrcr^ 
o 


Cj(0) * — (It /Bn) Rc 


- CD . 

\ cJ o*p (o /# *5 ») 


V ^ — fn / * v -- 


-- j* t -* t" i c r /••t/ % *• / • i 

* i \-i c l""4' r • ■*' * 


- exp [«,>,' (r . 8) 


H has been chosen such that a [wsltivu It means a divergent laser beam, and the In- 
tegration Is along the propagation p. th from s^) at th*> source to e « /. at t:.- recelvei 
From the Kolmogorov theory ot turbulence* 1 

^(o ; sj »8.10C N 2 0- ,,/3 ( 8 > 

Tlie gamma function is defined 


>(B. s). [2(4 -s)/ik] [(s-ika 2 )/(*-lka 2 >] 


whoro 


1 Ik 

r ? * “IT 


Combining equation 7 and equation 8 to obtain an expression for log-amplltudo variance 


c ,(°) - k 2 n«, 


jj d3C N 2 ^ a) | doo" n/0 ^exp |(o/2)ReJy (t, ajjj -c 
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I 


( 13 } 


For Or Q B 0 Y'-.. a) (f) 

C l* (0) r r <-5/6) ^RcJ | ds C N 2 (a) j[-l/2 Reyfz, s)] 5/ ~ [-l/2y (/., s) j 5/fl J j 


4 ~ r<- 5 / 6 ) k 7/r ’ Ho | j ds C N , 2 (s) )J 5/0 (1) 5/ G 


<H) 


( 15 ) 


For | 5/fl * 0. 259 + 0.905 i F (-ft/fl; « 0.08 

C, B CO) = O.«Ck 7/0 ( da C N 2 («) [(/•-») (ft 1 

* O 

The maximum vt»!ue of £(z- "> & occurs when s « z/2 indicating turbulence halfway 

between source anti receiver is most effective Ir. causing scIntlUmioa, when z »* 

2 

where s fij is the d* stance at which («)—•• O. 


rn 


Cj 8 (0) <* 0.56k 7/fl f m C n 2 (h) s 5/C ds 

J r 


<i«> 


ff the propagation is In a vertical direction to n receiver at altitude II and zenith angle O, 
the log-'impiftude varltirice becomes 

$ 1 ( g 

C N 2 (h) [(H-h> dh (1 7) 


Or for H » b where h is the altitude nt which 
m m 


V(h>- 


r h 

C^CO) -0.5(1 k 7/6 *oc 1I/6 0 \ C N 2 (h)h 5/6 dh 

* ft 


( 1 «> 
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It should be noted that there Is a reciprocity between plane wave space to earth 
propagation and spherical wave earth to space propagation. Because of this reciprocity, 
the scfntllhitlon of a star la the same as that of a ground-based laser source viewed 
from space, and equation 18 Is the same as Tatarskys' 1 equation 8.17. This fact was 
userl In the Cl DOS- If laser detector experiment In order to relate upward propagation 
to downward -propagation. Tills reciprocity Is also an Important factor tn the analysis 
of this experiment. 

lor noii|)oint sources (o^>0) the expression for log-amplitude variance Is more com- 
plex. Hi 1 3 problem has been analyzed by Fried** wlto finds that the sphcrlcnl-wave, Ing- 
nmplltude variance may be modified by n correcllon factor (Figure o — i) which Is a 
function of the normalized source size 0, 

2 

;i K k 6 o /h^sccfi (jgj 

The sources used In the experiment could not be made true point sources. Therefore, 
the normalized source slzo was made equal for both transmitters. This required that 
the standard deviations of beam profiles have the rnllo 


* 1.55 


where subscript a refers to the argon laser (X- o.f.i-IS x I0 _fi m) and *» to llie ciirlon- 
dloxlde laser (X B 10.6 x 10 ® m). The laser outputs were made C mm (argon) nr.d 
30 mm (carbon-dioxide) by pas.ifng them through optical systems which were adjusted 
to also provide a divergence of approximately one mllliradian. The optical transmitter 
system Is described in section 3.4. 1. 

Hie normalized source size for both transmitters for tho turbulence profllo of the ex- 
periment at a zenith angle of 45° was 


flO.fi x 

,o- C 

1 0. 51 45 

\ 10" 


fi “ «t£r 2 / h soce*^- 92 - * ltr> 

o o 1000 X 1. -I 


“ 0.381 


which resulted in an apo.turo averaging correction factor of 0.5. The balloon was flown 
at an altitude of 28.0 km above ground level, and was at a zenith angle of approximately 
<55° during the period data was taken. The transmitter was located on a level plain 1.2 
km above mean sea level. 




Using these factors and the turbulence profile of the previous section we get (altitude of 
the receiver was 28.0 km above ground level) the spherical wave log-amplitude varl- 


C, (o) = 0.56 k 7/C sec 11/6 O 
teO.56 k'^ G sec U/ ' 6 Q 




™ 2 . -1/3 -h/h .5/6 .. 
C.. h c o h ah 
No 


"1.10 k 7</6 [o. 89 C No 2 h o 3/2 J =< 1. 10 k 7/6 Jl. 18 X 10~ 9 J 

1.24 >: 19~ 9 x 7/0 


« 1.24 x 10 _9 x 1.81 x 10 8 « 0.224 (argon) 

= 1.24 x 10~ 9 x 5.44 x 10 6 = 0.00671 (carbon'dioxice) 
After correction for aperture averap ng 

Cj(0) » 0.59 (0,224) * 0.112 (argon) 

C L (0) * 0.50 (O.OOt'M) = 0.00337 (carbon-dioxide) 

To theso variances we must add th-5 scintillation caused by the tropopause 

5/6 

C.(0), ^O.GCk^w 11 ^ C + r<tl-hj h „ 1 


T 'trop 


P _JL 
II - 


- 1. 10 k" "* 4.3 x 10-“ [ <28 -° - 

- 6.62 X 10“ l ° k 7/C 

» 0.C2 x 10~ 10 x 1.81 x 10 8 o 0.120 (argon) (24a) 

» 6.62 x 10 -10 x 5 14 X 10 G *= 0.00360 (carbon-dioxide) (24b) 

The total log-amplitude variance duo to the turbulence near the ground plus the 
turbulence at the tropopause becomes 


Cj(0) = 0.112 + 0.123 *• 0.232 (argon) 

■ 0.00337 + 0.0036 ° 0.00697 (carbon-dioxide) 




» -W n • . 


& 


A pair of sensors aboard the balloon payload were used to monitor the lrradlance caused 
by the lasers. Instantaneous lrradlance levels were telemetered' to the ground and re- 
corded on magnetic tape. Later, the tape was played back In 30-second Intervals Into a 

2 

Hewlett-Packard correlator from which sample values of a were obtained. Normalised 

p 7 

power variances were converted to log-amplitude variances by equation 3. Hufnagel 
has predicted that the strength of scintillation as measured by log-amplitude variance 
should Itself be a log-normal parameter. Therefore, In C, (0) was plotted versus 
cumulative probability on normal graph paper to check the distribution of Cj (0). A 
normal distributee plotted (n this manner will appear as a straight line. Results are 
shown in Figures 0-5 and 6-6. The parameters of the distributions are listed below: 


0.5145 pm 

10,6/ ptr. 


In C,(o) -1.8 

-5.4 


<r Jn Cj (0) 0.6 

0.6 


The magnitude of C,(o) should, from equation 15, be proportional to the 7/6 power of the 

wavenumber. Therefore 



t C l (0) 3 0. fjI-15 pm = HO. 6 x 10' 6 I 
[CjfOlJ [6.5115 x 10- G J = 34.1 

10. G pm 


(26) 

From tho measured log-amplitude variances 



Median [(0,(0)] . .'J-* , 

Median [(0,(011 10>6#|m • “ 

36.9 

(27) 


The mean 


Cj (0) 

Cfi) 

C^oi 


from the data are* 






-1G.2 

c 

-5. 22 
o 


0.198 

0.0054 


(0.5145 pm) 
(10.6 pm) 


(28a) 

(28b) 


♦For a log-normal variable p = exp £k + Aj-J where p Is the mean of the log-normal 
variable and k and are the mean variance of the corresponding normal vnrlablea. 


* 



v * rr«0' 
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CUMULATING PROBABILITY - °£RCENT 

CUMULATIVE PROBABILITY* Of LOG-AMPLITUCE VARIANCE KOR C0 2 V.'AVELENGTH 
(10.6 mtcrc»T!«tcf») FOR A HOUR PERIOD. 4 30 - 7;lX) AM MST 21 ^CT08ER 1970 AT 
HOLLOMAN AFB, NEW MEXICO. BATE I EXPERIMENT 

DATA PLOTTED ON NORMAL PROBABILITY SCALE. STRAIGHT LINE INDICATES 
NORMAL DISTRIBUTION. 


Figure 6-6. Cumulative Probability of Log-amplitude Variance for CO,, Wavelength 


r ,. v< i4 

P- 


1 


Experimental values of Cj <0) determineo by optical measurement a^e compared to 
values computed from the turbulence profile In the following table: 


Cj(0) 0.5145 fim 
Cj(O) 10. C n m 


Experimental 

0.198 

0.00540 


Computed 

0.232 


0.00097 


Exp/Computed 


Hie central-limit theorem predicts that the distribution of irradlance should bo log- 
normal. l"hls was checked In Figures 6-7 and 6-8 by plotting the cumulative probability 
of the logarithm of a normalized signal on a normal distribution scale. The straight line 
behavior of Iwth of these graphs again indicates a log-normal distribution. 

In order to assure that the experimental values of variance were not causing a scatter 
of the results due to insufficient sample lengths, an analysis of the statistical error 
was performed. 

Standard statistical techniques provide r. means to estimate the rrfeet of finite data 
sample length. l or a number, N, of independent observations of a normally distributed 
random variable, x, it is known that the sample variance obeys a chi-square distri- 
bution with n degrees of freedom. 


JLf-1 = X 2 

2 A n 


n»N-l 

c 2 

ts sample variance 

2 

actual '\trianco 

2 

X R »* chi-square random variable 


of cid-squaro determine the 1 - a confidence Interval for sample variance, 


X *®/2 

!•# 


^<o 2 < 

X - 


X"; 1 -o/2 


(°$/S) Ml 
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CUMULATIVE PROBABILITY - PERCENT 

THE LOGARITHM OF NORMALIZED RECEIVED SIGNAl(S) FOR THE ARGON WAVELENGTH 
(5145 A) AS A FUNCTION OF CUMULATIVE PROBABILITY ON NORMAL PROBABILITY PAPER. 
t 49:05 MST 21 OCTOBER 1970 AT HOLLOMAN AFB, NEW MEXICO. BAPE I EXPERIMENT 

VARIANCE 0.72 

LOG-AMPLITUDE VARIANCE *»0.18 
So - MEAN SIGNAL LEVEL 


Figure 6-7. Logarithm of Normalized Received Signals for Argon Wavelength 
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CUMUlAirvE PROBABILITY - PERCENT 


THE LOGARITHM OF NORMALIZED RECEIVED SIGNAl(S) FOR THI: CO 2 WAVELENGTH 
00.6 mlcro»r«t«t) AS A FUNCTION OF CUMULATIVE PROBABILITY ON NORMAL 
PROBABILITY PAPER. 6:47:30 AM MST 21 OCTOBER 1 970 AT HOLLOMAN Af 8, NEW 
MEXICO. BAPE I EXPERIMENT 

VARIANCE C.04 

LOG-AMPLITUDE VARIANCE *0.01 
So - MEAN SIGNAL LEVEL 


Figure G-8. Logarithm of Normalized Received Signals for COg Wavelength 








that Is, thr actual variance o ^ Is known to fall within the above Interval with a confi- 
dence of 100 (1 - a) percent. The key parameter that determines the above Intervnl Is 
N, the number of independent samples. 

N*=2B c(j T (31) 

where \ 

B ~ equivalent white noise bandwidth 
eq 

T = record lengrh 


B In turn can oe estimated from 
eq 

B (°) 

n — 




2 / % ( y )d r 


where 


11 O') *• autocorrelation for random variable 


Topical values of 1!^ ,/j for nrRr>n and CO,, data lend to estimates for B and In turn N 
for a given sample record length, T. The value of chi-square for a particular N and is 
obtained from stardard tables. Results are given below for a 90-oerccpt confidence In- 
terval for each siipial ami a sample record length of 30 seconcs. 

Statistical Krror Analysis 


Data 

B eq < tJz > 

N 

a 



Argon 

150 

10,500 

0.1 

10,740 

10,203 

co 2 

35 

2,100 

0.1 

2,208 

1,995 


90 % confidence Intervals 
Argon 0. 93 S 2 ^ 2 ^ 1.C2 S 2 

C0 2 0. 95 S 2 <ffj 2 £ 1 . or. S 2 

Therefore the argon variance has error bars of * 2'e and CO 2 has * 5T for 90q confi- 
dence aial 30-second record length. 


• v? 


V*-.— “V f J \ 
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6.3 CORRKLAT i ON COE’FFICI ENT 
The covariance of Ir radiance is defined as 


c i (p) = «i 1 -T 1 >(I 2 -T 2 ) (1) 

v. hero {lie angle brackets denote an ensemnle average of ail |mIi n of (mints in the re- 
ceiving plane separated by a distance p. The variance o' irradlance Is simply tlr. 
variance for zero separation 

Cj (0) - 


In order to obtain a measure of the Independence of the scintillation nt two [mints In She 
receiving plane, the covnrtanc ia usually divided by the varla-'X to obtain the 
Ir radlanee-correlatli • coefficient 


y<fi ) 


C, <P) 

<7^ 


v 


<fi-rr> 


(2) 


Fit/. maur Ice' has shown, that the Irradlance-corrclatlnn confident may l>e related :o 
log-amplitude statistics 


y(p) 


e\p F-«C, <p|]-l 

<•:.[> f-iCj (0)J-1 


<.■** 


Whe.« Cj(p) •» <(1 J -Tj) (l 2 - !,» 

If ' .e define the log-amplitude correlation coefficient na F (p), 
. C. <p) 

then F{c , _ _t 

Cj (0) 

exp flCjiO) p<p)j -I 
V (P) “ exp'pe^Tj T 

For weak turbulence , 


(U 


y ( 0) 


From Fried , 




FW 


8. 16k 2 C N T(-|) 


it 


Re 


5/6 


(6) 


(7) 


(jW 7 


{.,8/3 r (Il/6)J 
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s 


m' 


l r -v5/0 r r ikf> 2 , , \ 

- 5 |i. (,-.»/*«> j « ,F,r-| ! i. r ^ 5 .]| 

n • 


( 8 ) 


For n propagation path over which C^, « 0 for all s except for a thin layer of turbulence 
of ihlcloiessAs << s this may be replaced by 

.,»vww r^- /3 , ,„Ui : 

C ! (0> ‘ “ 0e I ,«/:> 'JJ\ " 2 { k* f 


c 


2 F 

L \ 6 


1 1 [ •IfZ-Sim 


(o') 


Furthermore the turbulence need not bo constant nnd therefrre C^/'A a may bo replaced 

by 

C 2 

Cl, (s Ids 


The log-amplitude correlation coefficient- becomes 

3 /.I. 5/ <5 


( 10 ) 


F(f* « 


fcft y7 - s '"rn ■ -'<? il] " 2 [-ffia] m-ji » / V,H 1 “q£— ]| 


-f [-^j 


5/6 


He 


3.86 Ro|J -<-—■) 


5/3 


5/6 


Ref| 5/ti ) 2 5/3 f (1 1/G) 




( 11 ) 


which reduces to 


5/6 


FW - 3.8* «e{} -4.09 [-S^jyj] 


( 12 ) 


! F j *jp ,,f ~( T~ u) ~T ] 0 c <‘ nnu ®nl hy-pergeomclrtc 


The function 
function with an Imaginary argument. 


pr-srriiS '.f> 




' h r *aj ** \i ; s 


■rtfs® 
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t 


< 7 -»\ 7 . 
k s 


2 


JTl 


Then 'be hyporgcornotric function may he evaluated bv a Kunv era 3 function 


, F, (a; b; |x) 1 ♦ ♦ 


My (IV r 

(*>>„ 2 : 


(a) n (lx f 

(b) n n ! (H) 


where (n) « a (n * li (a * 2) (a * n - II ft*. « • 

** 3 

<b) » b (b ♦ 1) (b • 2| (b.o-I) (l» - I 

* o 

I t • •« Vk.j l la ut t lk.,« 

[-4 j -wo«> • T ;r* , -ji?r <« a 


i2J. *ur 


-JI22.— us . 

1 Mi. >-21 


r I -0. S33 t.N * 0.0.1I7X 2 . O.OiHSOIx 3 . O.tiJ.J*' 5 
- G.OO0O77ZX 1 ♦ . . . . 


Wo may normalise equation 12 by replacing p uy p / fkl.\ obtaining 


r (^=r) ' 3 - SG Ro {} - 4 - 09 f*) 




\y^* a 


I 


Using this equation, wc evaluate p / J\ L ' , »«{}. 

and \ as n function of x in the following:* 

Urv 


o 

Re j | 5 ' 6 F 

. -t . f v I ( 

F ( P 


Re ji jFj j^- T 

•'•'"JI 

(ynr 

a 

0.259 


1.000. 

0.500 

0.4 31 


0. 486 

0. 707 

0.603 


0.256 

1 . 000 

1 . 068 


0.020 

1.414 

1.870 


-0. 084 


The log-amplitude correlation corTicienl is sliown as n function of foe normalized 
parameter p/\/ XL' in Figure 6-3. The normalizing factor is known as the correlation 
length and corresponds to the separation fer which x = 1. 


Using the relationship of equation 6 botweeu y(p/ and F (° //Tp ) . wc caa oal- 

\ v : / V J / 

culate the theoretical dcrn.-ndei.ee of irradiance-corr.'lation coefficle t: 

y ('Val) 




o j 2 - 0. 20 

9 

Oj ~= 0 . 4 ; 

0 

0 

0 

0 

0.2 

0.2 

0.13 

0.1 

0 4 

0.4 

0.31 

0.22 

0 .G 

0.0 

0.50 

0.39 

0.8 

0.8 

3.73 

0.65 

1.0 

1.00 

1.00 

1.00 


In order to evaluate the actual correlation coefficient , we must analyze the turbulence 
as though It is made up of thin layers of turbulence of varying strengths and at varying 
distances from the source. 

Because the thermosonde measurement showed only one layer of turbulence and from 
other measurements we know that turbulence exists at the tropopause, we will use a 
two- layer model of the atmosphere, although the results of the analysis will make ap- 
parent how mor' complicated models can be handled. 

•See more complete tabio In Appendix B 
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For a single la; er of turbulence, 


Cj 8 (o) 


C n 2 As Re — 


2 t:/3 r (11/G) 


x 1 F 1 [ - / V* 1 '■ ] 

L M/.- s) z J 


(for s < < 7.) 




2 8/3 p (11/G) 


t |+j 


D = Ho j i 

l 


> F .[] 


c C . _ A s 


'N N 


Cj (p) - AC* 


■ [b.* 


* Bs 5/3 - cn* 5 ' 


c.V) = 'AC • CDs 5/<i 
1 N 


D * 0.2C> 


* - 0.2G A C* C b 
is 


C N ‘ '«*< * S/ ‘ 


- 0.26 A C N C 


For a two layer turbulence model , 


c l‘ <« " ac ni [ Bs , 5/ ’ - [ B "s 


>,<< » 2 


Dj = 0.26 


Tliereforo 


CjV) - a c N * |^n •*» - o. 26C« i ' ! ’ /6 J + ac n * 2 [-bs 2 5/3 - cd 2 s 2 5/6 j 


+ AC K2 I Bs 2 
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s 


S 5 *. 

Defining F(p) = Cj (p)/Cj (o), we have for the single layer model 


r “Vi • AC u [ ,ls i f/, ‘ ' c >1 


- 0.26 AC^ C 


CD 1 ~ ° 9 i 
0. 26 C 


5/6 


Fo~ ihe two layer model, 
F (P>l +2 


AC N1 0.2 6 c] + AC* 2 [bs 2 5/I1 - Ci)J 


0.26AC [C^ + C^] 


C N1 [ nB I ' " C D l] 


N2 


- ”-2« C [<Vl * C KsJ 


j»», 5/6 - C t>J (231 


L C N1 >> C N2 


F(P) 


1+2 


C N1 [ C D 1 - 

-, 5/r ’] * 

c + 

N2 

[c D 2 - Us., 5/S 

fc d - n s . 5/<5 

♦ 0.26 C C * 
M 

1 

^ + 
€„ n 

fc n - bs 5/0 ] 

L l 1 j 

I + 

N? 

L 2 ? J 

0.26 C 


c * 

^N1 

0.2G C 


F 4 C 


N2 


N1 


F (PL 


(24 ) 


U C N1 * C N2 


F <p) 


1+2 


C 


N1 


c + + c f 

C N1 ^N2 




F(P)j + 


c 

“ 

N2 


+ _ V 


C N1 4 C N2 



1 


F (P) 2 (25) 


» 




s 


and for weak scintillation 


1 

1-4 

* 55 
O 

V + 

Q + 

C N2 

1 c + + c + 

L N1 L N2 J 

y l 

C + *f c + 
C N1 ^N2 


' 1+2 

From this result, the method of analysis for multilayer models Is .pparcnt. 
k 


(26) 


F(P)= £ 
n * 1 


'N1 


L 


'NT 


J 


F(Pj) + 


C N2 


'NT 


F(P 2 ) + 


'Me 


'NT 


V 

. v 


F 0\) (271 


where C^ T Is Uie integrated value of cjj from 0 to * 

From Figure 6-10, the correlation length for a particular altitude can be obtained. 
Using this length, y^ and y^ can be calculated and, using the measured values of 
*"N1 al ’^ from the experiment, a theoretical shape for T (p) may be obtained. 
This ha« been -dune graphically In Figure 6-11 using the correlation coefficient curve 
of Figure C-2 and the following paragraphs 

CO 


Boundary Layer 
Tropopauso 


Mini. 

700 

’3,200 


Argon 

0.70 

0.12 


3.-3 

0.55 


c N y ^ L 

0.40 

0.60 


The experimental values of correlation coefficient are plotted In Figure 6-12 and re- 
corded In Table 6-2. Using the results of Figure 6-12 and converting F (p) to y 
through Table 6-2, the theoretical log-irrndlancc correlation coefficient for various 
values of log-amplitude variance are shown for comparison. 

The experimental and theoretical fit show excellent agreement indicating that the theory 
appears to give an accurate prediction of correlation coefficient. 


I 
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1 ' Tabic 6-2 

Tabic of Experimental Irradianco Correlation C 


Argon Data 

Carbon-DI 


Time (GMT) 

p(CM) 

y 

Time (GMT 

1 . 

13 29 16 

2 

0.949 

1 . 

13 32 40 

2. 

13 27 00 

2 

0. 939 

2. 

13 20 10 

3. 

.13 25 -15 

2 

0.921 

3. 

13 38 00 

4. 

13 23 00 

3 

0.811 

4. 

13 17 10 

5. 

i3 20 10 

4 

0.803 

5. 

13 08 30 

6. 

13 17 10 

6 

0.772 

6. 

13 46 15 

7. 

13 14 50 

8 

0.665 

7. 

13 48 45 

8. 

13 12 55 

12 

0.477 

8. 

13 09 30 

9. 

13 12 15 

14 

0.352 

9. 

1 ! 50 00 

to. 

13 08 30 

16 

0.400 

10. 

11 49 05 

12. 

[ 13 47 30 

22 

0. ?<2 

!!. 

11 P! 33 

12. 

13 09 30 

2C 

0.367 

12. 

11 56 50 

13. 

12 08 07 

35 

0.307 

! 13. 

13 53 20 





14. 

13 55 00 





lo. 

13 57 20 
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6.4 APKMTl’ItK A VERAGING OF SCINTILLATION 

Scintillation Is a p .oblcm mainly to laocr communications systems working In the 
visible or near infrared. Far-Infrared wavelengths, because of the 7/6 (xmc r of v/ave- 
mimlx*r dependence, are onliy slightly affected. There ..re only two methods to avoid 
sein'illation. The first, selection of a site which has low turbulence. Is partially 
effecti'T, but since nearly half of the scintillation Is due to tropopause which is too 
high to be affected by terrain changes. only small reductions can be made through site 
selection. Hulett 1 notes that astronomers, at many different locations a ound the 
world, observe nearly Identical levels of scintillation. 

The second method of reducing scInlUlntb'n is through the use of large apertures to 
average out the (n'enslty of the fluctuations. Since the ('tensity In ihc roceiv g plane, 
as has U-on shown in section 6.3. is not well correlated over separations larger ’hmi a 
few tens of cm, the aperture? a large receiver will l*e divided Into many small cells 
which scintillate independently from adjacent cella. Because of this behavior, the 
scintillation as measured by 


for aj.-'rturcB containing many independent cells should (mj proportional to ) A/n where 
n is the wiraUr of independent areas. Also, if on the average all of these areas arc the 
same sl/.e. the «/n* is proportional to the diameter of the receiver (I)) so ih.t for large 
apertures is proportional to 1/D. 

Figure 6-13 shows the ape Mure dependence of the average value of a as measured by 
Hufton^, Pro the roe"*, and Young*. For very small apertures at wavelengths In the blue 
anti green portion of the visible spectrum (A '- 0.500 f 0. 1 00 m) the magnitude of a ^ 

Is approximately 0.35. As the aperture Is Increased, decreases rapidly until 
an a|H*r»ure of at«out 25 cm (10 Inches! is reached. Beyond 25 cm, the rate of decrease 
in scintillation Is much slower, ami approaches the theoretical J/D dependence for 
large apertures. For D* 56 1 cm (20 Inches) 


r -r~ 


. , Vl ^>1' SX.,. . 
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6.5 AUTOCORRELATION OF IRRADrANCE AND POWER SPECTRAL DENSITY* 


The autocorrelation of Irradiance Is defined as 

R r (T) - <[I(t)-T] [l(t + T)-?]> = C r (r) (1) 


where I(t) is th~ irradiance at time t, i (t + r) is the irradiance at the s'mc point 
at time t + r, and T is the average irradiance. This is usually normalized by dl\‘.ding 
by the variance of Irradiance to obtain the irradiance-correlation coefficient 


>'( r ) 


C,(r) 

r. ■ — — 

C r (0) 


<fi(t)-r] [t(t* t)-!]> 

<u- 1)^> 


( 2 ) 


where the coefficient *s now a function of time delay T rather than of spatial dis- 
placement (p) as in equation 2 of audio.. C.3. For a single layer of turbulence at dis- 
tance s from the .source, a displacement p ' corresponds to a displacement 



(3) 


In the detector plane at ■/. from the projection shown in Figure 6- 14. If p is replaced 
by the correlation length in the detector plane, 


»■ -^[f] 


(<) 


■fho time necessary for a turbulon to traverse this distance is 



(5) 


where n is the wind velocity normal to the propagation path. 

Therefore If Taylor's hypothesis* of frozen-ln turbulence is assumed, T q is the time 

required for a turbulon In the plane of turbulence at s to traverse the correlation 

length («y5i£) in the receiver ptane, and is therefore simply the correlation time. 

Taylor's hypothosls is valid whenever v / XL<< L , that is when the correlation length 

0 2 

Is less than the outer scale of turbulence. From Fried , the outer scale of turbulence 
•Data processed only for the argon wavelength. 
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at 700 m altitude’ (lowest layer of turbulence) Is 40 m which is much larger than the 
correlation length for that altitude (0.70 m) and therefore the requirements of Taylor's 
hvpotehsls arc met for all the turbulence layers since correlation length decreases 
with altitude while outer scale increases. 

The wind speed corresponding to the tropopause (14.4 km) is 24.7 ms * and the corre- 
lation length corresponding to this height is 0.087 meters. Therefore from equation 5, 
the correlation time is 


T 

O 


0.087 

24.7 


r 14,4 x 1 1 

L 23.4 x 1 0-1 J 


= 2.17 x 10 


-3 


sec 


For the low altitude barrier layer turbulence, we get the altitude of the turbulence 
directly from thermosonde measurements and find ft to be approximately 700 meters, 
with a wind speed of 10.1 ms -1 (19. 0 knots). The correlation length for this altitude 
is approximately 0.70 meters, therefore 


7 

c 


0.70 f O . 7 0 xjn 3 "| 
10. i 1-23 . 4 x liP’ .1 


.-3 


In m;iklng this computation the velocity of the wind has been used without correction for 
Us azimuth since by fortuitous circumstance the wind vector for both layers of tur- 
bulence was precisely at right angles to the propagation path nt the time of the measure- 
ment (1313 GMT). Hadthewind been at another azimuth ( 0 W ) the normal component of 
the wind velocity could have been computed by the following relationship lrom Young 4 


sec 0- 


b 


*• tan 2 0, r 


sin 2 (0, 


- »*»] 


1/2 


whore 0 T Is the zenith angle of the target and Is the azdmuth angle. Since tho 
barrier layer turbulence and tropopause turbulcneo were of approximately the same 
strength, wc would expect tho correlation time to approximately 


T o " 


(2.08 2,17) x 10 


-3 


2. 13 x 10 sec. 


0-41 




j 

\ 

I 

I 

I 

) 

i 
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The actual correlation time from measurement (Figure 6-15) was 2.40 x 10 -3 seconds. 
The difference between predicted and measured correlation time is therefore 



(2.40 - 2.13) x 10 
2.4C x 10' 3 


11 . 2 % 


which is well wlti.in the accuracy of the wind speed measurement. 

The normalized power spectral density of' irradiances calculated from the auto- 
correlation of irradiance data by Fourier transform methods is shown in Figure G-16. 
Since the power spectral density of irradiance is proportinnul to wavenumber (k =^~), 
by scaling the frequency axis, the same curve may be used for other wavelengths. 
Carbon dioxide scintillation for example would have frequencies 

y 0 .5H. r > • 10~ ° 

10.6 x 10-C “ ' , “ £ 


times as great as the argon vavn.ength. The correlation curve can also b' used for 
carbon dioxide by multiplying tho time axis by 1/0.22. 
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Table C-3 

Log-Irradiauco Autocorrelation 


T 

y 

r 

y 

T 

1 v 

0 

1.00 

2.0 

0.17 

4.0 

- 0.04 

0.1 

0.99 

2.1 

0.14 

4.2 

- 0.03 

0.2 

0 . 9 S 

2.2 

0.10 

4.4 

- 0.02 

0.3 

0.97 

2.3 

0.07 

4.6 

0 

0.4 

0.94 

2.4 

0.04 

4.8 

0.02 

0.5 

0.91 

2.5 

0 . 02 

5.0 

0.02 

0 .G 

0.89 

2.6 

0 

5.2 

0.03 

0.7 

0 . 8.1 j 

2.7 

• 0. 02 ! 

i 

5.4 

0.03 

0.8 

0.80 

2.8 

- 0.03 

5.0 

0.03 

0.0 

0.75 

2.9 

- 0.05 

5.8 

0.03 

1.0 

0.69 

3.0 

- 0.06 

6.0 

0.03 

1.1 

0.65 

3.1 

- C .07 

6.2 

o . o :; 

1.2 

0.59 

3.2 

- 0 . U 7 

C .4 

0.03 

1.3 

: 0.54 

3.3 

- 0.07 

c.o 

0 . 0 .) 

1.4 

0.48 

3.4 

- t .. 03 

6.8 

1 

0.02 

1.5 

0.43 

3.5 

- 0.07 

7.0 

0.01 

1.0 

0.38 

3.0 

- 0.07 

7.2 

0 . O 1 

1.7 

0.32 

3.7 

- 0.06 

7.4 

0.01 

1.8 

0.28 

3.8 

- 0.03 

7.6 

0.01 

1.0 

0.23 

3.9 

- 9.05 

7.8 

0 


Wavelength- 0.5145 Time 13:12:55 T In roilUsoeoixla 
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7. CONCLUSIONS 

In tho previous sections of this report we have deacrttwrf live «»*[*• r I irw'nt and tnilunl 
It* result*. In this section we will attempt to synnpatre what ha* t*-en l«am«at, "nd 
what c.'tn l*“ concluded from the experiment results. 

f n Meet inn 6. I, Flgve fi-4, the vertical profile of torUdetwe (n the ftrat fl km of the 
ntmosplirra at the time of the experiment t# presented, together with a mathematical 
model useful for uiormlmallng the furtsdenre dUir|t«4tlon. \ke fled that ;Ho *<a!e 
height of the atmosphere la lowcr-ihan predicted b> Fried * try a factor of shout three to 


In section <1.2, the log-amplitude variance Is ana), ted, a«wt it t» hl«owi ,h*t the et length 
of nclnlllhitlnn can only !*• accounted tor by t *toi«g lurUdenct strive tl km On the 
tutals of (he wind profile , temperature profile ami correlation length, the layer ts pre- 
dicted to t« at the tropopnuse. The destruction of log-amplitude variance* t« *h»wi> to 
lie log-nonnal' with a standard deviation of O.fl, The w-vrlength dejx’rdence (k ' n | of 
Mcf-diliation ia cc.Ufrmed for the two wavelength* used in the (eat. the irradiaiwe dis- 
tribution is checked and found to t>c log-normal, 1 In- log-amplitude variance calcu- 
lated in the btsls of the profile of lurtmU-nce of section cl. I l* compared and fours) to 
Ire in agreement with the experimental values. 

In section fi.3. the log-amplitude cortelation coefflc lent t* taleolatrd on the nf 

the turlsiience model of section <5.1 and shown to t«e In agreement with the e .pci < mor 
taliy measured values. The correlation length fa found to lie proportional io i/l , ami 
a method of determining the correlation le.vgth for multilayer turUdence model* (• 
derived. 

Hascd ujwn tho result* ntiove, tho red-jctlon of aclnitltatloa using & Urgo *. vie scope (a 
calculated In section 6.4 and ccmpared to measurements by standard astronomical 
techniques. The dependence of MClntillQtloa frequency on wind speed la to* ted in tec- 
lion 0.5 and compared with tho measured data. Cortelation time nod power spectral 
deosiUtcs are shown to bo directly releied to wind *(>eod* in the various turbulent 
layer a. 

Tho analysis of tho data obtained from the experiment tends to confirm the theoretical 
predictions of scintillation based upon the turbulence- profile. Wo tiellcve that this (.* 
the first experiment (n which the optical effect* ot turbulence upon laser beams 
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propagating over '* verilct.! path have lieen tested In a controlled manner. We have 
found no i s*e iu «hu h the data contradicts the theoretical predict Iona, and have ol>- 
taine-it nM«.rv*,( (v good • ••nflrmatlonis of many of the assets of the theory. The exper- 
iment dwi h i* allowed in* to develop a techi ijue of measuring the xrrtlcal profile 
tt,r««i K ti a t-slt.w.n laxmehtsl ihcrmosonde ilex Ire. 

Our < xpe/).»i«nt ha* rseentlallv l***t«d S< hmelt/er‘s h.|ustlon ?.s for log-amplitude m- 
xai i.tm •&. Hetauei! the assumptions made for the log-amplitude covariance r’piatfona 
are largch ' t; - »*aw as those for phase covariance (f<fuit((nn ?.U) hr resuits of the 
♦pertinent 4end to 1 1 ltd c redence to the theory fot phase slellslii s *a well as fo - am- 
jdi' .•<«■ statistics- Ttrn reader sht-uld rex lets Hs hmel 1 1 *• r ’« (is |>er to see the tetnarl- alile 
edn.llarlt) c if M|in4uio T.s a.xd and dw fart that Ifiey are faith lutscd ufun Uie 

s.-.iue l« di derivation. 

home of the more imix.rtiint results of (he experiment Are shown In Ihe follow |ng table. 

T»>4«# i 

Nt>i|>*n of • i.ttw of Ihe Major IfesuUa of the HA!*K f r xpertmcnt 

1, V.-rtl. a! Profile of t’ * - l' ’ h' ,/3 e‘ h/ h o < 8 (h - h 1 C * 

furl silt m e | f lifure <i— 4) '* 0 

- 4 x lit rrs h ■- Jmto m 

■«i ti 

-■ ■4.3 x |o * ’* m* '* h " 13.2 x fo 3 nt 

vf» p 

2, Strength of rklfttliiallon 


a. Wavelength dependence (p. 6-14) 

.r 7/6 

ti. Dlslrltaiitnn (ITevrc-s fi-h olid d-9) 
c. Iair- A mplitude Variance (Mean) 

log-normal 

1) 0,5143 pm 'figure 6-6) 

0.17 

2) 10,6 pm < Figure 6-7) 

o.ooa 

d. Distribution of Sample Variance* 


1) 7ype ( figure# fl-0 and 6-7) 

Log-normal 

2) Standard Deviation (j*. 6-14) 

0.6 



t:-; 
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Table 7-1 

Synopsis of Some of the Major Hesults of Ihe HAPE 1 Exjierlment (Cent) 


3. Correlation of Intensity 

1/2 

k. Wavelength Dejiendcnc" of Correlation A 
Length (Figure 6-13) 

b. Correlation Length* 

1) 0.5H5 pm (Figure 6-13) 

2) 10.6 pm (Figure <1-1.1) 

c. Height Dependence of Correlation 
length (p, 6-26) 

•1. Aperture Averaging of Scintillation 

- 1 

a. Wavelength Dependence (p. 6-36) CCA 

b. Aperture Dependence (Lar-'o Apertures) CCD ' 
(Figure fi-H) 

c. Standard llevlatlon of Power (l-arg CCo * “f,~ 

A|'orturos) (p. 6-36) ^ 


22 Cm 
100 cm 



3, Correlation ‘l line 

a. Wind Sj«*etl Dei tend ence <i>. 6-44) 

b. Distance Dependence (p. 6-44) 

C» Wavelength Dependence (p. 6— M) 


-I 

0U- 

n 

CCh.'k 

CCA ,/2 


. J 


Let ua now consider how the result* of the experiment have anted laser communication 
theory. From the results, ve can estimate the lime dependence of po’ser received nt 
a ground receiver from a spacecraft, or the power received at a soacecraft from a 
ground receiver. \)'e can do this for any wavelength, for any sire receiver, for any 
zenith angle, awl (using rouulls from DAPF ID for any lime of day. Tills information 
can be directly applied to *110 estimation <•( bit error rnle. We can estimate 'he length 
of Revere drojwuta, their frequency, and the jiercent of time that n laser communica- 
tion* system will operate nt It* designed bit error rate or algnal-to-nolse ratio. He- 
cause of the multitude of factors affecting the system performance. It Is very <ltfflcult 
to present the result* In a single graph. Rather tlie results must be used to compute 
system performance for each new system a* It l« developed. Our new knowledge of the 
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atmospheric effects of turbulence should enable us to design laser systems which will 
operate successfully through atmospheric turbulence. 

One area which needs further investigation is the frequency and depth of severe inten- 
sity dropouts. These dropouts have little effect upon the means and variances but 
severely effect the bit error rate of pulsed binary laser communications systems. 
Because of the very high slgnal-to-noise ratio and dynamic range required to perform 
this type of mer 3urement, it was beyond the scope of this experimer;. 

A continuation of the BAPE I flight of 1970 was performed in the fall oi 197J . This ex- 
periment, BAPE II, has added to the data of this report (particularly in the area of 
turbulence profiles) and a report on its results is under preparation. 

Tha BAPE system is the only existing equipment capable of measuring the statistics of 
the optical effects of turbulence upon laser beams propagating from earth to epace. 
While the BAPE I and II flights have provided a great deni of tniormation, many inter- 
esting and important experiments remain to be do-e. The most Important ot these are 
the direct measurement of the phase statistics of optical propagation, and the measure- 
ment of the frequency of severe amplitude dropouts. Work in these areas should be the 
next problems of laser propagation to be investigated. 
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APPENDIX B 

EVALUATION OF LOG-AMPLITUDE CORRELATION COEFFICIENT 
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